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The continuous pumping of atoms into a Bose-Einstein 
condensate via spontaneous emission from a thermal reservoir 
is analyzed. We consider the case of atoms with a three-level 
A scheme, in which one of the atomic transitions has a very 
much shorter life-time than the other one. We found that 
in such scenario the photon reabsorption in dense clouds can 
be considered negligible. If in addition inelastic processes can 
be neglected, we find that optical pumping can be used to 
continuously load and refill Bose-Einstein condensates, i.e. 
provides a possible way to achieve a continuous atom laser. 



I. INTRODUCTION 

Duringthe last years, the fruitful combination of laser- 
cooling M and evaporative cooling has allowed the 
experimental achievement of the Bose-Einstein conden- 
sation (BEC) in trapped weakly interacting alkali gases. 
Such remarkable achievement has stimulated an enor- 
mous interest || . The subsequent efforts mainly concen- 
trated on two different areas. On one side, the BEC offers 
an extraordinary opportunity to test condensed matter 
and low-temperature phenomena. In this respect, very 
recently several striking results have been reported con- 
cerning superfluidity phenomena |^,^| and generation and 
dynamics of vortices HQ and dark-solitons ||[)| . On the 
other hand, the macroscopically occupied matter wave 
can be manipulated by atom optical elements, that can 
be combined to provide new tools for precision experi- 
ments . Besides passive optical elements recently also 
active elements, that provide phase coherent gain have 
been demonstrated |0. Also a new field, called Non- 
Linear Atom Optics (NLAO), has rapidly developed dur- 
ing the last years. Several remarkable experiments have 
been recently reported in this reflection of BEC 

from an optical mirror [fl2[ and four-wave mixing ||l3| of 
matter waves. 

Among the results related to BEC and NLAO, one of 
the most important experimental achievements was the 
realization of an Atom Laser. As a coherent source of 
matter waves, the atom laser will lead to new applica- 
tions in atom optics. Its impact in the field is compa- 
rable to the one of light lasers in light optics. The first 
realization of an atom laser was achieved via pulsed rf- 
outcoupling from a BEC The coherent character of 
the source was demonstrated in a landmark experiment 
jlql) in which two atom-laser pulses where overlapped, 



showing a clear interference pattern. 

Since this first realization, several groups have build 
atom lasers using (quasi-) continuous outcoupling from 
the BEC, either by using rf fields 16 lj], or by employing 
Raman pulses [Q. However, the continuous outcoupling 
represents just a half way towards a cw atom laser. The 
continuous loading of the condensate still remains to be 
incorporated in experiments. Without a continuous re- 
filling of the BEC, the atom laser output lasts only as 
long as some atoms in the BEC are kept. Like in the 
development of light lasers the availability of cw atom 
lasers would open the way to " high power" and precision 
applications. 

Two different physical mechanisms could provide a 
continuous pumping into a condensate. On one hand, the 
collisional mechanisms , in which two non-condensed 
atoms from a reservoir collide, and as a result one is 
pumped into the condensate, whereas the other carries 
most of the energy and is evaporated. On the other 
hand, the optical pumping of reservoir atoms into a BEC 
via spontaneous emission processes has been also pro- 
posed ||(J. If this reservoir could be filled in a (quasi-) 
continuous way by laser cooling techniques, one would 
benefit from the large cooling efficiency of laser cooling 
compared to evaporative cooling, allowing for a consider- 
able increase in atomic flux produced by an atom laser. 
For the latter, it is crucial that the spontaneously emit- 
ted photons cannot be reabsorbed, because otherwise a 
heating is introduced in the system, and BEC can be 
neither achieved nor maintained pl[ |. 

Several possible dynamical and geometrical solutions 
for the reabsorption problem have been proposed dur- 
ing the last few years. The geometrical proposals are 
based on the reduction of the dimensionality of the traps 
|^3|,|lJ. It is easy to understand that assuming that the 
reabsorption cross section for trapped atoms is the same 
as in free space, i.e. ~ l/k 2 L , the significance of reabsorp- 
tions increases with the dimensionality, in such a way 
that the reabsorptions should not cause any problem in 
one dimension, have to be carefully considered in two di- 
mensions, and forbid condensation in three dimensions. 
Therefore, cigar-shape and disc-shape traps have been 
suggested. However even severe deformations of the trap 
do not allow more than modest reductions of the reab- 
sorption heating j24]. Other suggestion consists in us- 
ing a strongly confining trap with a frequency u> ~ ujr. 
In this case, it has been proved |2f| that in two atom 
systems the relative role of reabsorption in such a trap 
can be significantly reduced. It is, however, not clear 
whether this result would hold for many atom systems. 
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Another promising remedy against reabsorption heating 
employs the dependence of the reabsorption probability 
for trapped atoms on the fluorescence rate 7, which can 
be adjusted at will in dark state cooling In particu- 
lar, in the interesting regime in which 7 is much smaller 
than the trap frequency u>, i.e. in the so called Festina 
Lente limit pq| , the reabsorption processes, in which the 
atoms change energy and undergo heating, are practi- 
cally completely suppressed. However due to the slow 
time constants in this approach the cooling efficiency is 
greatly reduced. Another reabsorption remedy could be 
provided by the destructive quantum interference of the 
negative effects of the photon reabsorption in the, so- 
called, Bosonic- Accumulation Regime [£7| . 

In this paper, we concentrate on the continuous optical 
pumping into a BEC. The paper is divided in two differ- 
ent parts. In the first part we present a new scenario 
in which the reabsorption is suppressed, in much less re- 
strictive conditions than that of Festina Lente regime, i.e. 
without reducing the cooling efficiency. In this scenario, 
an atom possesses an accessible three level A scheme, in 
which one of the atomic transitions decays much faster 
than the other. By employing Master Equation (ME) 
techniques, we show that the dangerous reabsorptions of 
photons on the slow transition are suppressed because 
the respective coherences are destroyed due to the de- 
cay via the fast transition. By dangerous we mean here 
those reabsorptions that may lead to undesired heating 
of the system. Since the dangerous reabsorptions are not 
present, this scheme can be employed to continuously 
pump atoms into the lower level of the slower transition. 
In the second part of the paper we analyze the dynamics 
of such pumping in the presence of atom-atom collisions, 
and show that the combination of elastic collisions (evap- 
orative cooling), and bosonic enhancement of the spon- 
taneous emission, can create a condensate, and refill it 
in the presence of outcoupling or losses. Therefore, this 
scheme could be considered as a possible way towards 
a continuously loaded atom laser. As a possible exper- 
imental realization we consider laser cooled Chromium, 
but the ideas can be generalized to any atom that pro- 
vides an asymmetric three-level system. 

The structure of the paper is as follows. In Sec. |l| 
we introduce the physical model, as well as the quantum 
ME that determines the loading dynamics. In Sec. |l| 
we introduce the so-called Branching Ratio Expansion 
(BRE), which allows us to analyze the hierarchy of pro- 
cesses which occur in the system. In Sec. IV we analyze 
in detail the suppression of the reabsorption effects. Sec. 
[v| is devoted to the treatment of the atom-atom colli- 
sions. In Sec. VI we present the numerical Monte Carlo 
results of the loading dy nam ics. Finally, we summarize 
some conclusions in Sec. VII. 



II. MODEL 

We consider the case of atoms with an accessible three 
level A system (see Fig. I), formed by the levels |r), |e) 
and \g). The atoms are trapped in an isotropic harmonic 
trap which, depending of the internal state of the atoms, 
has frequencies uj ri uo e and uj gi respectively. This could 
be, for instance, the case of Chromium 52 Cr, in which 
the electronic levels would be 7 S3, 7 P4, and 5 D4, respec- 
tively. The transition |r) <-> |e) is assumed to be driven 
by a laser, which has a Rabi frequency Q. The spon- 
taneous emission frequencies associated with the tran- 
sitions |r) <-> |e) and \g) <-> |e) are, respectively, "f er 
and 7e 9 , such that 7 er 3> 7 eff . In the case of 52 Cr, 
7 er = 2tt x 5MHz > 7 eg = 2tt x 30Hz. The branch- 
ing ratio e = jeg/jer is therefore very small (~ I0~ J in 
52 Cr). The fact that e <c 1, will lead to the suppression 
of the reabsorption of scattered photons. 




FIG. f . Atomic scheme considered throughout the paper. 

In this section, we shall not consider the collisions be- 
tween the atoms in the \g) state. Such collisions are in- 
troduced in our formalism in Sec. [v|. In the following 
we take h = c = 1 for simplicity. Let us introduce the 
annihilation and creation operators of atoms in the |r), 
\e) and \g) states and in the trap levels s, I, and m which 
we shall call r s ,r|, ei,e\, and g m ,9m- These operators 
fulfill the standard bosonic commutation relations. 

The Hamiltonian which describes the coupling of the 
system of bosons to the laser field, as well as to the vac- 
uum electromagnetic modes is of the form: 



H = H + Her + H% + K) + Hf, 

with the following terms: 



(f) 



• Free atomic Hamiltonian (describing internal and 
center-of-mass degrees of freedom): 



(2) 



with wj, tof and lu^, denoting the energies of the 
level s of the \r) trap, the level I of the |e) trap, 
and the level m of the \g) trap, respectively, loq is 
the transition frequency between \r) and |e). 
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• Interactions of the laser quasiresonant with the 
transition |r) <-> |e): 



n 



,e- iWLt elr s +h.c. 



(3) 



where tj; iS is the Franck-Condon factor which de- 
scribes the transition between a level s of the |r) 
trap, and a level / of the |e) trap, and uul is the 
frequency of the applied laser. 

• Spontaneous emission processes |e) — > |r): 

£«/=-*ee/^ 



2 £ o(27r)3 (del "^) 



xms(k)e]r s a^ + h.c. 



(4) 



where where a,- , at are the annihilation and the 
creation operators of a vacuum photon character- 
ized by a wavevector k and a polarization /i, with 
polarization vector eg ; d er is the dipole vector of 
the transition |r) |e). 

• Spontaneous emission |e) — > |.g): 



^/ 'EE/'^ ^ 'A,,' 

XVlmWe^ma^ + /i.e., (5) 



where d e9 is the dipole vector of the transition 

IS) - 

• Free Hamiltonian of the electromagnetic (EM) field 
Hf = E/ (6) 



P = -iH e ffp + ipH e ff + Jp, 



(7) 



where p is the density matrix, H e f f = H^ff + H^X and 
J = J (0) + J (1 \ with 

F S = E W ^ + EK - 6 + 

s I 

E W ™#m5m +H er , (8a) 
H eff = ~^ e 9 E E a lmm'('4wl' e l'i ( 8b ) 

jWp = 2 ler J2J2^n>rUi'p4rs, (8c) 
= 2 7eff ^ a^^g^eypejgrn. (8d) 

/,m I' ,rn' 

Here i5 = — is the detuning, 



(9) 



whereas a^^;/ = a^ s , ; , + ia/ as ,j„ /3 imm ^ = $? mm , v + 

iPLm'l" with 

= J dkwCk)ms(k e 9 k)vv m ,(k e 9 k), (10a) 
-1„ f°° „ £ 3 



"lmm'1' = P / 



dfcW(fc) 



(10b) 



In the previous equations k indicates the solid angle co- 
ordinates, W(fc) represents the dipole pattern, k^ 9 is the 
wave number associated with the transition \g) <-> |e), 
and P denotes the Cauchy Principal part. Finally, /3^ sll , 
has an identical form as o^ im ///, but for the transition 
\r) <-> |e). 



Starting from the Hamiltonian ([j]), one can trace the 
full density matrix of the system over the vacuum modes 
of the EM field. Using standard methods of quantum 
stochastic processes Eqj, one can derive then the Quan- 
tum Master Equation in Born-Markov approximation 
p9j. In principle such ME fully describes all the pro- 
cesses which happen in the system, including eventual 
reabsorptions in the fast transition e) — > \r). However, 
for simplicity of the analysis we shall consider the case 
in which we can neglect the reabsorption phenomena for 
the reservoir atoms (|e) and |r)). Although what fol- 
lows is true also in presence of those reabsorptions, the 
approximation will allow us to concentrate on the much 
simpler problem of the pumping of a single atom from 
the reservoir {|e), \r}} into the \g) trap, where of course, 
collective phenomena are important, and therefore taken 
into account. In this case, the ME takes the form 



III. BRANCHING RATIO EXPANSION 

In this section we shall show that if e <C 1, the re- 
absorption effects on the slow transition can be safely 
neglected, because they occur with a probability 1/e 
times smaller than the spontaneous emission processes 
l e ) ~~ * Is) without any reabsorption. In order to demon- 
strate that, we shall perform the analysis of the different 
processes that could happen in the system. 

We shall consider the situation in which some atoms 
are already accumulated in the level \g) and the time to a 
single atom is being pumped from the level \r) to the level 
|e). This atom undergoes then the spontaneous emission 
process, and its further fate may be twofold. First, it may 
undergo transition |e) — > |r), and the emitted photon 
will then leave the system (we assume no reabsorptions 
at |e) «-> |r) line due to the low density of \r) atoms). 
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Second, the excited atom may undergo a transition |e) — > 
\g). We allow for the possibility that the emitted photon 
in this process may be reabsorbed several times, until 
it leaves the system or until the emission e) — > \r) will 
take place. We shall analyze the hierarchy of possible 
processes which can be produced on the time scale of 
the above described processes. We shall assume that no 
other atom is pumped from the level |r) to |e) within this 
time scale. Technically, this approximation means that 
we exclude the possibility of multiple quantum jumps 
from |e) to |r), i.e. we assume that the whole process 
consists of the sequence of processes, each involving an 
atom being pumped from |r) to |e), which then undergoes 
spontaneous emission processes (including reabsorptions) 
until it either lands on the level \r) or \g), after which 
another atom is being pumped from |r) to |e), and so on. 

This approximation is performed here just for reasons 
of technical simplicity, but we want to stress that: 

• The approximation describes well the experimental 
situation with Chromium atoms in which the Rabi 
frequency £1 < j er . The pumping process has thus 
indeed an incoherent character, consisting in a se- 
quence of jumps |e) — > \r) followed by spontaneous 
emission events. Of course, it may happen that sev- 
eral atoms are being excited simultaneously to |e). 
Each of them, however, behaves independently of 
the others, so that the analysis pertaining to just 
one excited atom is valid. 

• The result obtained below is indeed more general, 
because at any time scale the hierarchy of proba- 
bilities is maintained. The latter statement means 
that our results hold also for > 7 er . This can 
be at best understood using a dressed-state pic- 
ture with respect to the |e) <-> |r) transition. In 
the limit £1 3> 7 er the dynamics reduces to a sit- 
uation in which a single atom is being pumped to 
one of the dressed states |+) (|— )), and undergoes 
spontaneous emission consisting in arbitrary num- 
ber of incoherent |+) — > |+) (|— ) — > |— )) transi- 
tions, ending either in \g), or in |— ) (|+)). This 
process is followed by an arrival of another atom in 
the state |+) (|— )), followed by successful sponta- 
neous emission, etc. Each one of these steps can be 
understood with the model presented below. 

The formal solution of the ME (fjj) after a photon es- 
capes from the system is given by: 

/>oo 

Poo= / dtJie-^'f'tp}. (11) 
Jo 

We have used in the expression ( |il"| ) the shortened nota- 
tion£ e //p= —iHeffp+ipHlff. In the following we shall 

denote C^/ f p = ~iH { ^ f p + ipH^f f , C { ^ f p = -iB [ ^ s p + 

ipH^ff Since j eg -C j e r,Q, we are going to perform 
an expansion in the branching ratio e. We consider an 



initial state of the system given by |?/>o) = IV'o) ® IV'o' 7 ')) 
with |i/>o) = | -Wo j N±, . . ., N m , . . .), with N m denoting the 
initial population of the m-th state of the \g) trap. IV'o''") 
denotes the initial state of the |r) and |e) traps. We 
are interested in the probability to obtain a final state 

l^y> = l^> 

The branching ratio expansion must take into account 
the bosonic enhancement effects. Due to the bosonic en- 
hancement, and the fact that we consider that the ground 
state of the \g) can be macroscopically populated, the ex- 
pansion must be done in the parameter eiVo, and not only 
in e. Nevertheless, we expect that the expansion and the 
conclusions that we draw from it will remain valid un- 
til eNo — 1/2. To use this expansion in the considered 
experimental situation we have to limit ourselves to the 
case of iVo < 10 5 /2. That means, however, that the 
expansion can nevertheless be safely used during the on- 
set of the condensation, and the dangerous reabsorbtions 
can be neglected in that regime. The situation in which 
cNq > 1 requires the treatment of higher order terms in 
the expansion, but does not mean that in such a case 
the reabsorbtions will cause troubles. In the case when 
very many atoms are already condensed, the dynamics 
is dominated by bosonic statistical effects. The use of 
similar ideas and techniques as those employed in the 
Boson Accumulation Regime (BAR) ^71 should be then 
possible. 

Employing standard methods of time-dependent per- 
turbation theory in the small parameter eNo we obtain: 

(iMpoolV/) = A o + A 1 + A 2 + O((eN ) 3 ), (12) 

where Aj is the term of order 0((eiVo) J '). Let us now 
analyze step by step each of the terms of the branching 
ratio expansion (BRE): 

A. Zeroth Order 

The zeroth order term is of the form: 

r°° (o) 
A = <V/| / dtJ^ie^fftpoWf), (13) 
Jo 

This process is by far the most probable one, and im- 
plies an spontaneous emission into the \r) state These 
processes (and the subsequent repumping by the laser) 
induce a thermal distribution of the \e) and |r) traps, 
and do not affect directly the populations of the \g) trap. 

B. First Order 

The first order term of the BRE is of the form A\ = 
Aia + Ai b , where 

Aia = (ipf \ ! dtJ^ie^ff'poWf), (14) 
Jo 
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dt'e 



.(0) 
-e// 



(t-t 



''{^{e^/po}}}^/}. (15) 



Ai a corresponds to the case in which an atom in the 
state |e) decays (without any further reabsorption) into 
a state of the \g) trap (see Fig. ||(a)). An, is given by 
the interference of two different process: a process like 
the one considered in the zeroth order, and a process in 
which (i) a decay is produced into some state of the \g) 
trap, (ii) a subsequent reabsorption is produced from the 
same state of the |<?) trap, and finally a process as that 
of the zeroth order is produced. The processes described 
by Au, although containing reabsorptions, do not change 
the population distribution of the \g) trap, and can be 
considered as small quantitative corrections to the zeroth 
order term. 
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FIG. 2. Scheme of different processes which appear in the 
BRE. The solid line implies decays or reabsorptions, whereas 
the dashed line represents the photon which is spontaneously 
emitted and further reabsorbed, (a) represents the process 
of first BRE order A\ a ; (b) and (c) represent respectively the 
negative and positive reabsorption effects which appear in the 
second BRE order processes A2 a \ the interference between the 
processes (a), and (d,e) constitutes the term A?,b of the second 
order of the BRE. 



C. Second Order 

The second order term of the BRE is of the form A2 = 
A 2a + A 2b , with 



dt'e c °ff (t_t ' ) 



A 2a = (ip f \ I dtjM 

dt"e^'-^{c%{e^' PQ ))\\) (16a) 

A 2b = {ipf\ / dtj^ 
Jo 



(16b) 



Let us consider the term A 2a . This term involves pro- 
cesses in which an atom originally in some state of the e) 
trap, decays into some state m of the \g) trap, producing 
an spontaneously emitted photon, which is reabsorbed 
by other atom in some other state n^mof the lg) trap. 
These processes are of order e 2 , except the case in which 
n = or m = (Figs. ||(b) and (c)); in such a case, if 
the system is already condensed, the probability associ- 
ated with these processes is of order e 2 7Vo. We must note 
that the process of Fig. ^(b) introduces a negative effect 
of the photon reabsorption in our system, because pro- 
duces a non-condensed atom, while destroys an already 
condensed one (of course the opposite process of Fig. ||(c) 
corresponds to positive effects, and is of the same order). 
The term A 2b is due to the interference effects between 
the process considered in A\ a and the processes of Figs. 
||(d) and (e). These process do not cause any negative or 
positive effects of the reabsorption, and simply introduce 
small quantitative corrections to A\ a . 

As observed, the "bad" reabsorption processes which 
change the \g) trap population distribution (i.e. may lead 
to undesired heating) are of order e 2 No, and therefore are 
1/e times less probable than the single spontaneous emis- 
sion into the \g) trap without any reabsorption. Hence, 
the reabsorption effects can be safely neglected, i.e. the 
atoms in the \g) trap can be considered as transparent for 
the spontaneously emitted photons on the |e) — > |<7) tran- 
sition. We shall show in Sec. VI that this effect can be 
used to optically pump atoms from the reservoir {|r}, |e)} 
into the \g) trap, and eventually into a condensate cre- 
ated in it. 



IV. SUPPRESSION OF THE REABSORPTION 
EFFECTS 

In this section we analyze in detail the physical effect 
behind the decreasing of the reabsorption probability in 
the considered system. The physical picture can be un- 
derstood by taking a closer look to the expression (16a) 
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which, after eliminating the terms which do not change 
the populations of the \g) trap, takes the form 



V ft f-t' 

dt I dt' / di"(V/|i7 (0) 



{m-Am-})\ + w-xwwm iv>/>> 



where 



I//. ,/\\ (t-i') zj(l) -iffi?- t'i , v 



(17) 



(18) 



Therefore, the process can be divided into three parts: 
(i) From time to some t' , the system evolves following 

the effective Hamiltonian -£fg//> ( u ) ^ a spontaneous 
emission occurs from |e) to \g), followed by a reabsorp- 
tion; (iii) The system undergoes after t' the same dynam- 
ics as in the part (i), until time t where a jump is 
produced into the \r) state. 

As observed in the expression (|l7|), the term Ai a de- 
pends on the correlation of the amplitudes of probability 
of two processes (i— iii) in which (ii) is produced at two 
different times, t' and t" < t'. In the interval t" to t' 
a jump into r) is produced with large probability, and 
therefore the probability amplitude is reduced, roughly 
speaking, by a factor exp(— -f er (t' — t")). Therefore, the 
correlation decays very rapidly (in comparison with ^y~ g ) 
with the time difference t' — t" . This leads to the strong 
reduction of the reabsorption probability. 

A very intuitive picture of the underlying physics in 
3-level atoms can be obtained in the following way: The 
reabsorption on the \g) to |e) transition has an oscilla- 
tor strength j eg . However the excited state coherence 
is rapidly destroyed due to the decay of state |e) to |r), 
with a rate 7 er . Therefore the absorption line width on 
the \g) to |e) transition is dominated by j er . Similar to 
other broadening mechanisms the effective reabsorption 
cross section is reduced by a factor e = jer/jeg- The 
atomic sample can be e-fold more dense than in the nor- 
mal two level case before it becomes optically thick and 
reabsorption becomes a relevant process. 

It is perhaps interesting to mention here the differences 
and similarities with the Festina Lente regime [^6| , which 
constitutes a known way to avoid the reabsorption prob- 
lem. As pointed out previously, the reabsorption prob- 
ability is determined (as already shown in p6| ) by the 
correlation at different times of the amplitude of reab- 
sorption probability. In the case considered in p6| , the 
correlation decays with the same frequency T as the spon- 
taneous emission frequency of the system, and therefore 
the probability of decay plus reabsorption is that of the 
decay without further reabsorption (multiplied by some 
geometrical factor, specially important in the so-called 
free-space limit |E6[). In the case of the Festina Lente 
regime (r <C oj) fl26f the reabsorptions which do not pre- 
serve the energy are suppressed due to a different reason 
than that considered in the present paper. In the Festina 
Lente conditions the interference terms in the amplitude 
of probability of "bad" reabsorptions at different times 



have a phase which rapidly oscillates in time; therefore, 
the time integration leads to a strong reduction of the re- 
absorption probability (by the small factor T/oj). In this 
sense, therefore, the Festina Lente Regime can be under- 
stood in terms of an expansion similar to the BRE. The 
case considered in this paper is, however, different. The 
spontaneous emission rate "f eg is now a factor e smaller 
than the rate of decay of the correlation, given by 7 er . 
This explains why the reabsorption is a factor e less prob- 
able than a decay without reabsorption. 

As a final remark, let us point out that the BRE does 
not necessarly imply small spontaneous emission rates 
(7 eg < uj), as it was the case of Festina Lente, but just 
the relative condition j eg <C 7er- In particular, j eg could 
be larger (or even much larger) than the trap frequency. 
Therefore, in principle, the atoms could be pumped into 
the \g) trap much faster than in the Festina Lente limit. 
This is of special interest when considering a sufficiently 
fast continuous loading of a condensate. 



V. COLLISIONS 

In this section we introduce the collisions between 
atoms in the \g) trap. In particular, no collision is consid- 
ered between the atoms in the \g) trap, and those in the 
|e) and \r) trap. Such approximation is valid assuming 
the situation in which the reservoir atoms are at much 
larger temperature and lower densities than those atoms 
in the \g) trap. Since we shall consider a \g) trap of a 
finite depth, the eventual collisions with the relatively 
much hotter reservoir atoms would lead to losses in the 
\g) trap, and eventually to a depopulation of the conden- 
sate created in it. Since we consider a loss mechanism 



from the condensate anyway (Sec. VI) such collisional 
losses could be easily taken into account phenomenolog- 
ically in our simulations as an effective outcoupling rate. 

The effects of the elastic collisions within the \g) trap 
are accounted by a new term in the Hamiltonian ([!]): 



H, 



coll 



mi ,777 2 ,7773,7774 



^m 1 ,rn2,m 3 ,rn i gln i gln 3 9m 2 9rn 1 - (19) 



In the regime we want to study, only the s-wave scatter- 
ing is important, and then: 



U, 



4irh 2 a 



mi ,7712 )T7l3 ,7774 



R3 



d 3 x Ct^^m., (20) 



where ipm denotes the harmonic oscillator wavefunctions 
and a denotes the scattering length. 

In the following we are going to work in the so-called 
weak-condensation regime, where no mean-field effects 
are considered. This means that we consider that the 
mean-field energy provided by the collisions is smaller 
than the oscillator energy. Therefore the model is only 
valid to describe either the onset of the condensation, or 
the full condensation but with the constraint that the 



G 



condensate cannot contain more than, say 500-1000 par- 
ticles. A realistic calculation of the dynamics beyond 
the weak-condensation regime would require the self- 
consistent calculation of the atomic states, which due 
to mean-field effects would change during the dynamics. 
Such calculation could be possible by using Bogolyubov- 
de Gennes formalism, and will be the subject of further 
investigation. In this paper we concentrate in the weak- 
condensation regime, where one can apply the formal- 
ism of Quantum-Boltzmann Master Equation (QBME) 
pfj| , pT| to treat the collisional effects. By using similar 
arguments as those employed in the context of collective 
laser cooling in the presence of atomic collisions in the 
weak-condensation regime [^2|, one can show that the 
ME which describes in this regime the loading dynamics 
from the reservoir {|e), |r)} to the \g) can be divided into 
two independent parts: 



pit) = C coU p(t) + Ci oad p{t) 



(21) 



where Ci oa dp{t) is the rhs of Eq. ([7]), and C co u describes 
the collisions, and has the form of a QBME as that de- 
scribed in Refs. 

Summarizing, the dynamics of the system splits into 
two parts, (i) collisional part, described by a QBME, 
and (ii) loading part, described by the same ME (Q) as 
without collisions. The independence of both dynam- 
ics, constitutes the main technical advantage of consid- 
ering the weak-condensation regime, and allows for an 
easy simulation of the loading process in the presence of 
Is)— Is) collisions. In particular, we simulate both dynam- 
ics using Monte Carlo methods, combining the numerical 
method of Ref. (3lj , with simulations similar as those al- 
ready presented in Refs. |29f| . 

The numerical simulation of the collisional dynamics 
for a three-dimensional harmonic trap is demanding, due 
to both the degeneracy of the levels, and the difficulties to 
obtain reliable values for the integrals U[n\, n-2, 77,3, 77,4). 
Therefore, we shall limit ourselves to the use of the er- 
godic approximation |^l| , p3[ , i.e. we shall assume that 
states with the same energy are equally populated. The 
populations of the degenerate energy levels equalize on a 
time scale much faster than the collisions between levels 
of different energies, and than the loading typical time. 
Following Ref. |33|] the probability of a collision of two 
atoms in energy shells mi and 1112, to give two atoms 
in shells 7713 and 7774 (where this collision is assumed to 
change the energy distribution function), is of the form: 

P(mi,TO2 — > 777,3,7714) = A(mj + l)(rrij + 2) x 



N mi (N v 



Qrrii ~\~ ^m3,m 4 ) 



SmiSm2Sm 3 Sni4 

where g mk = (m^ + l)(wfc + 2)/2 is the degeneracy of 
the energy shell rrik, rrij — min{mi, 7772, 7773, 7774}, an d 
A = (Aa 2 uj g m)/ (ttTi). In the following we shall use for 
the calculations the mass of 52 Cr, and, since a is not 
known for 52 Cr, we shall adopt a scattering length a = 6 



nm (similar to that of Rubidium). We shall consider a 
trap frequency lo„ — 2tt x 1kHz. 



VI. NUMERICAL RESULTS 

In this section we consider two different physical prob- 
lems. First, we shall show that under appropriate con- 
ditions it is possible to load an initially empty trap for 
atoms in the state |s) via spontaneous emission from a 
thermal reservoir of particles e). This allows to achieve 
the condensation in the trap |s) in a finite time, which 
depends on the physical parameters. Secondly, we shall 
show that in presence of an outcoupling (or, as discussed 
in Sec. [v], in presence of trap losses) it is possible to 
maintain the condensate population by refilling the con- 
densate via the spontaneous emission from the thermal 
cloud. 

The zero order term of the BRE corresponds to the 
fastest transition, but does not affect directly the loading 
process. Here we shall concentrate on the process pro- 
vided by Ai a , i- e - the spontaneous emission |e) — * \g). 
The loading rate of a state \g,n) from a state \e,l) is 
therefore given by 

T(l,n) = 2 leg J dkW(k)\(n\e rkf \l)\ 2 (N n + l). (23) 

Let us assume that N ex atoms are in the level |e) with a 
thermal distribution of temperature T. The total loading 
rate into the state |<7, 77,) is provided by: 



r(n) = 5>(0r(i,n), 



(24) 



where p(l) = N ex exp{—fku e l/kBT}/ Z is the thermal 
distribution, with I = l x + l y + l z and Z = (1 — 
exp{— hu) e /kBT})~ 3 is the canonical partition function. 
Thus, the loading rate can be rewritten as: 



r(n) = 2 7eg / dkW(k) 



(77 e 



5> 



,-huiJ/k B T 



e-^ r \n)(N n + l) (25) 



Let us consider now the following conditions which 
greatly simplify the numerical calculation of the loading 
process. We are going to assume that u> e < uj g + 2aj rec /3, 
so that there is always at least a state of the excited- 
state trap in an interval ±cj rec (with u rec being the recoil 
/^■frequency) around any considered state of the ground- 
state trap. The ground state trap has a finite depth, 
i.e. it possesses a maximal energy shell at haj g m max 
p4j| . Therefore, only those |e) atoms with energies 
E < tuv g m max + Tiio rec = E max can effectively load the 
trap. Let us consider that the reservoir has a tempera- 
ture T such that ksT E max Q. Under the previous 
conditions, the expression between the brackets in (|25|) 



7 



can be safely substituted by 1/ ' Z. Since the previous con- 
ditions also imply hu e <C ksT, the loading rate can be 
rewritten as: 



r(n)= 7e// (7V„ + l) 



(26) 



with j e ff = 2^ eg N ex (huj e /kBT) 3 . In other words, if 
the previous conditions are satisfied, all the levels of the 
ground state trap are equally loaded. 

In our numerical simulations we consider a \g) trap 
which is cut at an energy shell m max = 10-60 (which for 
ui g = 2ir x 1kHz, implies a trap depth of 0.5-3.0/iK). 
In the simulations we have used a "virtual" trap of 
m-max + 10 energy shells in order to take into account 
the atoms which do not decay into the trap, and also 
to calculate the evaporative cooling dynamics. We have 
checked that such chosen "virtual" trap does not affect 
the physics of the problem. The levels m > m max are 
continuously emptied in the simulations, and those atoms 
are considered as lost. 



40000 60000 

t(l/co ) 




FIG. 3. (a) Condensate fraction (No/N) as a function of 
time, in units of the , for the case of j e ff = 0.01w s , 
a — 6nm and m max = 50; (b) Dynamics of No for the same 
situation; (c) Evolution of the energy per particle for the same 
conditions. 



We have numerically simulated the loading dynam- 
ics of an initially empty \g) trap for different values of 
the parameter Jeff, rn max and the scattering length. 
It is worthy to note that Jeff is related to the phase 
space density of the thermal cloud by the relation j e f / = 
2j eg x 5.2n ex \{T) 3 , being A(T) = (2nh 2 /Mk B Tf/ 2 the 
thermal de Broglie wavelength, and n ex the density in the 
|e) trap. Since the |e) atoms are considered at a tempera- 
ture far above the critical temperature of the onset of the 
condensation, then n ex \(T) 3 is necessarly much smaller 
than 1, restricting the possible range of values that j e ff 
can take. We study below the dependence of the loading 
dynamics on j e ff- It is also interesting to point out the 
limits of the large-temperature approximation, in which 
all the levels are equally loaded. If such approximation 
is valid, for a fixed phase space density <f>, a constraint 
to the density of the reservoir is introduced by n ex S> 
(j)(MLorh max /2Trh) 3/2 , with m max = m max +L0 r e C /uj g . For 
the case of 52 Cr and uj g = 2ir x 1kHz, this constraint im- 
plies n > 7A54>mUax x 10 n cm~ 3 . 

As an example of loading dynamics we show in Fig. 
|3| the case of j e ff = 0.01cj 5 , m max = 50, and uj g — 
2tt x 1kHz. For the case of 52 Cr, 2j eg = 200s" 1 , this 
would imply a phase space density (f> — 10~ 5 , and a 
constraint n ex 10 8 cm~ 3 in order to satisfy the large- 
temperature approximation. Due to the random charac- 
ter of the process, we have averaged over several Monte 
Carlo realizations, in order to obtain smoother curves. 
Typically the loading process is characterized by a time 
scale of the onset of condensation, after which, as ob- 
served in Fig. ||, the condensation in the ground-state 
of the trap appears. In the case of Fig. || the onset of 
the condensation appears at approximately 3 x lO 4 ^;" 1 , 
which for oj g = 2ir x 1kHz would require t = 4.8s. After 
a time 10 5 u> g (16s) more than 1000 atoms are condensed 
with a condensate fraction of 70%. 

Let us discuss the physics behind the presented nu- 
merical results. Since the trap is initially empty, and all 
the loading rates are equal, the higher shells of the trap 
are initially loaded with larger probability, due to their 
larger degeneracy. Once a sufficient number of particles 
is loaded inside the trap, the collisions produce the evap- 
oration of part of the particles, whereas part of them are 
transported via collisions to lower levels of the trap. Dur- 
ing this process the energy per particle is continuously de- 
creased in the trap (see Fig. |^(c)). Finally, the condensa- 
tion is produced in the ground state. Once this happens, 
a new mechanism which reduces the energy per particle in 
the trap appears, namely the bosonic enhancement of the 
spontaneous emission into the condensate. In effect, the 
condensate loading becomes faster. When the number of 
condensed particles becomes large the QBME equation 
is no more valid, as discussed in Sec. [V|. The mean field 
effects appearing beyond the weak condensation should 
not, however, distort the qualitative effect of speeding 
up of the loading rate . It is also worthy to note that 
for condensate densities larger than 10 15 atoms/cm 3 in- 
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elastic processes (as three-body recombination) are ex- 
pected. Such processes would lead to losses of condensed 
atoms, which can be repaired by the continuous loading 
in a similar way as described below for the case of atom 
laser outcoupling. Let us also note that, eventually, if the 
number of condensed particles were much larger than the 
number of levels in the trap, the system would enter into 
the so-called Bosonic Accumulation Regime, where the 
vast majority of the decays were produced into the con- 
densed state. 

With this physical picture in mind, it is possible to 
understand the dependences on the different physical pa- 
rameters. From the previous discussion it becomes clear 
the crucial role played by the collisions in the process. 
Therefore the larger the collisional probability, the faster 
the condensation is achieved. This can be obtained in 
a two-fold way, (i) by increasing the scattering length a 
and/or (ii) by increasing the number of trapped atoms. 
Point (i) have been studied in Fig. ||[ where we have 
analyzed the case of m max = 50, j e ff = 0.01w g for val- 
ues of the scattering length ranging from a = 1.25nm to 
a = 24nm. As pointed out previously, larger values of 
the scattering length produce faster condensation onset. 
Note, however, that the time of onset of the condensa- 
tion reaches a constant value for large scattering length; 
this is simply because there is, as pointed above, a time 
in which the initially empty trap is being loaded, and 
in this initial time the collisions play no significant role. 
Such initial time depends basically on m max and 7 e t / , as 
pointed out below. 

A faster increase of the number of particles in the trap 
can be achieved in two different ways: 

• Increasing 7 e //- In Fig. |5| we show the example of 
Timax = 50, and a scattering length a = 6nm, for 
j e ff/ujg ranging between 1 and 0.01. As one could 
expect the onset of the condensation appears faster 
for faster loading. In order to achieve a larger 7 e //, 
it is necessary either to increase the phase space 
density of the thermal cloud, or perhaps more in- 
terestingly to consider experimental situations with 
bigger j eg . Note that an increasement of -f eg by a 
factor 100 in the case consider above, still very well 
satisfies the requirements of the BRE. We should 
stress once more at this point that one of the main 
advantages of the BRE is the fact that contrary 
to Festina Lente it is possible to pump with a fre- 
quency larger than that of the loaded trap, with- 
out having problems with the photon reabsorption. 
Therefore, it could be in principle possible to con- 
sider faster natural transitions, or even controllable 
quenched ones. In such fast onset of the 
condensation could be possible for low phase space 
densities of the reservoir. 

• Increasing m max . In this way the total number 
of atoms in the trap becomes larger at any given 
time. As an example, we have considered in Fig. 
H the case of j e ff = 0.01o7 g , for m max = 30 and 



60. The increasement of the condensate fraction 
becomes faster for larger m max , although the onset 
time remains approximately the same (Fig. ||(a)). 
The absolute number of condensed atoms is larger 
for larger values of m max (Fig. ||(b)), because the 
total number of atoms is also larger. Therefore 
for practical purposes (e.g. detection) it would be 
recommendable to use larger m max - However, we 
must stress at this point that a larger m max needs 
a larger reservoir density n ex if one wants to remain 
within the limits of the large-temperature approx- 
imation. 




5 10 15 20 25 



scattering length (nm) 

FIG. 4. Time of onset of the condensation for 
j e ff = 0.0lLo g , rrimax = 50, as a function of the scattering 
length. 
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FIG. 5. Time of onset of the condensation for a — 6nm, 
m m ax = 50, as a function of 7e//- 
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20000 40000 60000 80000 1e+05 

t(l /to,) 

FIG. 6. (a) Condensate fraction (Nq/N) as a function of 
time, in units of uij , for the case of Jeff — 0.01cj 9 , and 
a = 6nm, and m m ax = 60 (solid) and 30 (dashed); (b) Dy- 
namics of iVo for the same situation; (c) Evolution of the 
energy per particle for the same conditions. 



We have also analyzed the case in which an already 
formed condensate is emptied via outcoupling, and con- 
tinuously pumped via spontaneous emission from the 
thermal reservoir. In our simulations we just consider 
outcoupling from the condensate, although similar meth- 
ods could be employed to simulate losses affecting the 
whole trap. We simulate without outcoupling the cre- 
ation of a condensate as described above, for the case of 
7 e // = 6.28o; g , and m max = 10. For the case of 52 Cr, 
and 2j eg = 200s~ 1 . This represents a quite large phase 
space density 6 x 10~ 3 ; however, we must again stress 
that the BRE allows to work with much larger j eg , and 
therefore with much lower phase space densities of the 
reservoir. At t = 350ms (when Nq ~ 950), we begin 



the outcoupling. We have analyzed different outcoupling 
rates j out (Fig. 0) , and monitored the population N af- 
ter 16s. This allows us to find a critical threshold £o (for 
this case 1.14) for the ratio £ = jout/jeff- For £ > £ the 
loading ("gain") is faster than the outcoupling ("loss"), 
and the number of condensate atoms increases with time. 
For £ < £o) the number of particles decreases and stabi- 
lizes for a lower N . For £ <C £o no condensate can be 
kept. 




FIG. 7. Number of condensed particles after 16s (after ini- 
tially loading an empty \g) trap during 350ms) for the case 



of lu — 2ty x 1kHz, m ri 
outcoupling rates j out 



10, 



7e// 



6.28(^0, and different 
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FIG. 

the case of m max 
Jout = (1.17 



Condensate population, averaged during 40s for 



6.28w 



2tt x 1kHz, 



10, J eff 

/(t))7e//, with < f(t) < 0.05 chosen ran- 
domly from an uniform distribution. 



It turns out to be important to maintain the popula- 
tion as constant as possible, for the reasons that we clar- 
ify below, and therefore to work in the regime of £ = £o- 
It is however not an easy task, due to the stochastic na- 
ture of both, the collisions and the pumping mechanism. 
In order to stabilize the noise optimally, i.e. to preserve 
the population of the condensate as constant as possi- 
ble it is useful to introduce a random temporal variation 
of the outcoupling rate. Fig. || shows the averaged dis- 
tribution of population of the condensate during 40 s of 
continuous outcoupling, for the case of an outcoupling 
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rate lout = (1.17 - /(*))7e»//, with < /(t) < 0.05 cho- 
sen randomly from an uniform distribution, and the same 
conditions as in Fig. [?]. The population of the conden- 
sate is maintained quasiconstant with an average value 
of (N ) = 940, and a variance ({N§) - {N ) 2 ) 1/2 = 80. 
During these 40s, 3 x 10 5 atoms are extracted from the 
condensate, with a rate of 7500 atoms/s. 

Let us briefly comment about the importance of keep- 
ing the population of the condensate as constant as pos- 
sible. The mean-field interaction translates the varia- 
tions of the condensate population into variations of the 
energy of the outcoupled atoms, being the variance of 
the energy related to the variance of the condensate den- 
sity: cr{E) = 4:TTh 2 aa(no)/m. Therefore, the narrower 
the population distribution of the condensate, the more 
"monochromatic" will be the atom-laser source, and con- 



sequently the larger the coherence time will be 36 1. Let 
us point out finally, that an additional way to control the 
fluctuations of the condensate population could be pro- 
vided by monitoring the energy of the outcoupled atoms, 
which would inform about the variations of the conden- 
sate density. Such information could be used in a feed- 
back loop to dynamically adapt the outcoupling rate to 
reduce the energy variance of the outcoupled atoms, and 
therefore increase their temporal coherence. 



VII. CONCLUSIONS 

In this paper we have analyzed a possible mechanism 
which could allow the creation and continuous loading of 
a condensate from a thermal reservoir, by optical pump- 
ing. In order to achieve such loading mechanism, it is nec- 
essary to guarantee that the reabsorptions of the sponta- 
neously emitted photons do not lead to undesired heating 
of the atoms in the trap. We have analyzed a particular 
scheme which allows to satisfy such condition. In this 
scheme an atom forms a three level A system, in which 
one of the transitions decays much faster than the other 
one. By using quantum Master Equation techniques we 
have shown that the very small branching ratio between 
both transitions induces very large reduction of prob- 
ability of the reabsorption processes which change the 
population in the lowest state of the slower transition. 
We have explained such effect by identifying the pho- 
ton reabsorption as a process whose probability depends 
on the correlation between the reabsorption amplitudes 
at different times. Such correlation is rapidly destroyed 
by the fast decay into the other possible channel. The 
destruction of this correlation causes the desired effect, 
i.e. the reduction of the "bad" reabsorption processes, 
responsible for possible heating. 

Once we have shown that the reabsorption has no sig- 
nificative effect on the system, we have analyzed the 
loading dynamics from a thermal reservoir, using Monte 
Carlo simulations, including the atom-atom collisions in 
the QBME formalism. We have analyzed the loading of 



an initially empty trap, demonstrating that the onset of 
the condensation appears after a finite time, which de- 
pends on the physical parameters of the system. The 
condensation appears due to the joint combination of 
thermalization via collisions, evaporative cooling due to 
the finite depth of the considered trap, and bosonic en- 
hancement of the pumping process. We have also an- 
alyzed the continuous refilling of the condensate, once 
it has been formed, taking at the same time into ac- 
count continuously outcoupling. We have shown that 
the refilling mechanism allows the compensation of the 
losses introduced by the outcoupling, and we have ana- 
lyzed the best strategies to keep the condensate popula- 
tion quasiconstant, which is important in order to achieve 
a "monochromatic" atom laser output. In the paper we 
have only analyzed the outcoupling mechanism, but the 
same reasonings applies to possible condensate losses, 
produced by inelastic processes, such as three-body re- 
combination, or collisions with the thermal atoms in the 
reservoir. 

All our simulations and estimates have been done for 
Chromium atoms, and for the parameters of the experi- 
ment currently performed at the University of Stuttgart. 
It is however interesting to stress that the same scheme is 
general, and in particular can be applied for other atomic 
systems, such as Magnesium [p7| . As a final remark, 
we would like to stress that the mechanism of avoiding 
the "bad" reabsorption processes, considered in this pa- 
per (i.e. regime of BRE) allows for faster pumping than 
other reabsorption remedies (such as Festina Lente, for 
instance) , and therefore allows for more effective compen- 
sation of the condensate losses. It offers a novel and in- 
teresting perspective towards a continuously loaded atom 
laser. 

We acknowledge support from Deutsche Forschungs- 
gemeinschaft (SFB 407), from the EU through the TMR 
network ERBXTCT96-0002, and from ESF PESC Pro- 
gramm BEC2000+. We acknowledge fruitful discussions 
with W. Ertmer, E. Rasel, K. Sengstock, and M. Wilkens. 



[1] S. Chu, Nobel Lecture, Rev. Mod. Phys. 70, 685 (1998), 

C. Cohen-Tannoudji, Nobel Lecture, ibid., 707; W. D. 
Phillips, Nobel Lecture, ibid., 721. 

[2] M. H. Anderson et at, Science 269, 198 (1995); K.B. 
Davis, M. O. Mewes, M. R. Andrews, N. J. van Drutten, 

D. S. Durfee, D. M. Kurn and W. Ketterle, Phys. Rev. 
Lett. 75, 3969 (1995); C. C. Bradley, C. A. Sackett, and 
R. G. Hulet, ibid. 78, 985 (1997). 

[3] Proceedings of the International School of Physics "En- 
rico Fermi", Course CXL, eds. M. Inguscio, S. Stringari, 
and C. Wieman (IOS Press, Amsterdam 1999). 

[4] C. Raman, M. Kohl, R. Onofrio, D. S. Durfee, C. E. 
Kuklewicz, Z. Hadzibabic, and W. Ketterle, Phy. Rev. 



11 



Lett. 83, 2502 (1999). 

0. M. Morago, S. A. Hopkins, J. Arlt, E. Hodby, G. [26] 
Hechenblaikner, and C. J. Foot, Phys. Rev. Lett. 84, 

2056 (2000). [27] 
M. R. Matthews, B. P. Anderson, P. C. Haljan, D. S. 
Hall, C. E. Wieman, and E. A. Cornell, Phys. Rev. Lett. [28] 
83, 2498 (1999). 

K. W. Madison, F. Chevy, W. Wohlleben, and J. Dal- 
ibard, Phys. Rev. Lett. 84, 806 (2000). 
J. Denschlag, J. E. Simsarian, D. L. Feder, Charles W. 
Clark, L. A. Collins, J. Cubizolles, L. Deng, E. W. Ha- [29] 
gley, K. Helmerson, W. P. Reinhardt, S. L. Rolston, B. 

1. Schneider, and W. D. Phillips, Science 287, 97 (2000). [30] 
S. Burger, K. Bongs, S. Dettmer, W. Ertmer, K. Sen- 
gstock, A. Sanpera, G. V. Shlyapnikov, and M. Lewen- [31] 
stein, Phys. Rev. Lett. 83, 5198 (1999). 

C. S. Adams, M. Sigel, J. Mlynek, Phys. Rep. 240, 145 [32] 
(1994). 

S. Inouye, T. Pfau, S. Gupta, A. P. Chikkatur, A. Grlitz, [33] 

D. E. Pritchard, W. Ketterle, Nature 402, 641 (1999). 

M. Kozuma, Y. Suzuki, Y. Torii, T. Sugiura, T. Kuga, [34] 

E. W. Hagley, L. Deng, Science 286, 2309 (1999). 

K. Bongs, S. Burger, G. Birkl, K. Sengstock, W. Ertmer, 
K. Rzazewski, A. Sanpera, M. Lewenstein, Phys. Rev. 
Lett. 83, 3577 (1999). 

L. Deng, E. W. Hagley, J. wen, M. Trippenbach, Y. Band, 

P. S. Julienne, J. E. Simsarian, K. Helmerson, S. L. Ros- 

ton, and W. D. Phillips, Nature 398, 218 (1999). 

M.-O. Mewes, M. R. Andrews, D. M. Kurn, D. S. Durfee, [35] 

C. G. Townsend, and W. Ketterle, Phys. Rev. Lett. 78, 

582 (1997). 

M. R. Andrews, C. G. Townsend, H.-J. Miesner, D. S. 
Durfee, D. M. Kurn, and W. Ketterle, Science 275, 637 [36] 
(1997). 

I. Bloch, T. W. Hansen, and T. Esslinger, Nature 403, [37] 
166 (2000). 

S. Burger, K. Bongs, K. Sengstock and W. Ertmer, Proc. 
Int. Scholl of Quant. Elec, 27th Course, Erice, Italy 
(1999), in print Kluver Academic, Amsterdam, 2000. 
E. W. Hagley, L. Deng, M. Kozuma, J. Wen, K. Helmer- 
son, S. L. Rolston, and W. D. Phillips, Science 283, 1706- 
1709 (1999). 

for a recent proposal E. Mandonnet, A. Minguzzi, R. 
Dum, I. Carusotto, Y. Castin, J. Dalibard , Euro. Phys. 
J. D (to be published); |cond-mat /9909378 and references 
therein. 

R. J. C. Spreeuw, T. Pfau, U. Janicke and M. Wilkens, 
Europhys. Lett. 32, 469 (1995); U. Janicke and M. 
Wilkens, Adv. At. Mol. Opt. Phys. 41, 261 (1999) and 
references therein . 

M. 01shan'ii,Y. Castin, and J. Dalibard, Proc. 12 th Int. 
Conf. on Laser Spectroscopy, M. Inguscio, M. Allegrini, 
and A. Lasso, Eds. (World Scientific, Singapour, 1996). 
A. M. Smith and K. Burnett, J. Opt. Soc. Am. B 9, 1256 
(1992). 

T. Pfau, and J. Mlynek, OSA Trends in Optics and Pho- 
tonics Series 7, 33 (1997). 

Y. Castin, J. I. Cirac, and M. Lewenstein, Phys. Rev. 
Lett. 80, 5305 (1998). 

U. Janicke and M. Wilkens, Europhys. Lett. 35, 561 



(1996). 

J. I. Cirac, M. Lewenstein, and P. Zoller, Europhys. Lett. 
35, 647 (1996). 

J. I. Cirac and M. Lewenstein, Phys. Rev. A 53, 2466 
(1996). 

C. Gardiner, Handbook of Stochastic Methods, Springer 
Verlag, 1985; C. W. Gardiner and P.Zoller, Quantum 
Noise, Springer- Verlag, 1999; Howard Carmichael, An 
Open Systems Approach to Quantum Optics, Springer- 
Verlag, 1993. 

L. Santos and M. Lewenstein, Phys. Rev. A, 60, 3851 
(1999). 

C. W. Gardiner and P. Zoller, Phys, Rev. A 55, 2902 
(1997) 

D. Jaksch, C. W. Gardiner and P. Zoller, Phys. Rev. A 
56, 575 (1997). 

L. Santos and M. Lewenstein, Appl. Phys. B 69, 363 
(1999). 

M. Holland, J. Williams and J. Cooper, Phys. Rev. A 55, 
3670 (1997). 

The large-temperature approximation is assumed since 
it simplifies the numerical simulation. Although this ap- 
proximation has certain physical consequences, we want 
to stress that similar effects as the ones discussed in this 
paper are not only possible for the case of lower temper- 
atures, but should in principle allow for faster loading. In 
such a case, the Franck-Condon factors of the different 
transitions should be properly calculated (c.f. po|p9| ). 
For a discussion in the case of Festina Lente regime see 
L. Santos, Z. Idziaszek, J. I. Cirac, and M. Lewenstein, 
in print in the Speccial Issue of J. Phys. B on "Coherent 
Matter Waves", ed. K. Burnett (2000). 
M. Holland, K. Burnett, C. Gardiner, J. I. Cirac, and P. 
Zoller, Phys. Rev. A 54, R1757 (1996). 
The latter possibility is considered by the W. Ertmer/E. 
Rasel group at the University of Hannover. 



12 



